Quantum effects after decoherence in a quenched phase transition 
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We study a quantum mechanical toy model that mimics some features of a quenched phase 
transition. Both by virtue of a time-dependent Hamiltonian or by changing the temperature 
of the bath we are able to show that even after classicalization has been reached, the sys- 
tem may display quantum behaviour again. We explain this behaviour in terms of simple 
non-linear analysis and estimate relevant time scales that match the results of numerical sim- 
ulations of the master-equation. This opens new possibilities both in the study of quantum 
effects in non-equilibrium phase transitions and in general time-dependent problems where 
quantum effects may be relevant even after decoherence has been completed. 
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I. INTRODUCTION 

The emergence of classical behaviour in quantum sys- 
tems is a topic of great interest for both conceptual and 
experimental reasons [1]. It is well established by now 
that the interaction between a quantum system and an 
external environment can lead to its classicalization; de- 
coherence and the occurrence of classical correlations 
being the main features of this process (for a recent 
overview see [2]). 

A seemingly unrelated physical problem where the in- 
teraction between a main system and its surrounding en- 
vironment is central is in determining the dynamics of 
a phase transition. Usually, a change in the properties 
of the system or the bath, forces the system to change 
phase via an out-equilibrium evolution. It is natural to 
ask what role decoherence plays in the phase transition 
and conversely, how the time dependent nature of the 
process affects the classicalization of the system. 

In this article we explore two concurrent avenues. We 
look at what may happen with the decoherence process 
when we have a time dependent setting (so far this prob- 
lem has been mostly studied in kicked or driven systems; 
see for example Refs. [3] and [4]). This is a very gen- 
eral question, and we use to guide us a simple toy model 
that naturally includes time dependent features. This 
model also happens to mimic some properties of a non- 
equilibrium second-order phase transition, giving us some 
clues as to what may happen in a realistic case. 

The paper is organised as follows. In the next section 
we introduce our model and review the physical role of 
the different terms in the relevant evolution equations. 
We describe how the 'phase transition' is implemented 
and discuss estimates for the different time scales in- 
volved. In Section III we present the results of a series 



of numerical simulations for the evolution of an initial 
configuration of two de-localised Gaussian wave packets. 
This system is subject to a sudden quench via an instan- 
taneous change in the frequency sign. Both the cases 
where the temperature of the environment is kept fixed 
and allowed to change at the quench time are studied. 
We support the numerical results with a detailed ana- 
lytical analysis. Section IV contains similar results this 
time taking as initial condition a single Gaussian state 
centred at the global minimum. In Section V we discuss 
the time dependent evolution of the linear entropy in the 
model, illustrating the loss of purity of the system and 
clarifying the physical nature of the results previously ob- 
tained. Section VI contains final remarks and the main 
conclusions of the paper. 



II. THE MODEL 

We will start by considering a quantum anharmonic 
oscillator coupled to an environment composed of an infi- 
nite set of harmonic oscillators. The total classical action 
for the system is given by: 



S[x. q n ] = S[x] + S[q n ] + S int [x, q„ 



ds 



C n xq n , (1) 



where x and q n are the coordinates of the particle and 
the oscillators respectively. The quantum anharmonic os- 
cillator is coupled linearly to each oscillator in the bath 
with strength C n . This coupling leads to a simple quan- 
tum Brownian motion (QBM) model commonly used in 
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the study of the quantum to classical transition [5,6]. 
Tracing over the degrees of freedom of the environment 
one obtains a master equation for the reduced density 
matrix of the system. From this one can derive the fol- 
lowing evolution equation for the corresponding Wigner 
function [2]: 

W t (x,p,t) = {H syst , W r } PB - ^xd 3 ppp W r 
+ 2 7 (t)d p ( P w T ) + D{t)dl p w, - f(t)d 2 px W r , (2) 



where 
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7 (i) is the dissipation coefficient, D(t) and f(t) are the 
diffusion coefficients. r)(t) and v(t), the dissipation and 
noise kernels, are given respectively by: 



77(f) = / du>I((j) sin uit 
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where I(lo) is the spectral density of the environment. 

The first term on the right-hand side of Eq.(2) is the 
Poisson bracket, corresponding to the usual classical evo- 
lution. The second term includes the quantum correction 
(we have set h = 1). The last three terms describe dissi- 
pation and diffusion effects due to coupling to the envi- 
ronment. In order to simplify the problem, we consider 
a high-temperature ohmic (I(u>) ~ u>) environment. In 
this approximation the coefficients in Eq.(2) become con- 
stants: 7(f) = 70, / ~ 1/T, and D = 2j a k B T. The nor- 
mal diffusion coefficient D is the term responsible for de- 
cohercncc effects and at high temperatures is much larger 
than 7q and /. Therefore in Eq.(2), we may neglect the 
dissipation and the anomalous diffusion terms against the 
normal diffusion. It is important to note that the high- 
temperature approximation is well defined only after a 
time-scale of the order of l/Ck^T) ~ 70/-D (with h = I). 
The relevant period of evolution for our systems takes 
place at times comfortably larger than this time-scale, 
safely in the validity regime of the approximation. 

Time dependence will be introduced in the Hamilto- 
nian by imposing a sudden change of sign of fig (typical 
quench) . This mass term is taken to be positive initially, 
the original symmetry being broken by f2 2 becoming neg- 
ative. On a second stage we will also consider the case 
where the temperature of the environment T changes 
with time. The change in the potential leads to the for- 
mation of degenerate minima mimicking the breaking of 



symmetry in a second order-phase transition. In a realis- 
tic model one should address this problem in the context 
of quantum field theory [7] . This is an extremely difficult 
problem since non-perturbative and non-Gaussian effects 
are relevant in the dynamical evolution of the order pa- 
rameter undergoing the transition and clearly numerical 
simulations are out of the question. We trust that any 
non-trivial type of behaviour that may be a feature of our 
simple quantum mechanical model will also be present 
(and likely more strongly so) in the infinite dimensional 
case. 



III. DE-LOCALISED INITIAL STATES: 
QUANTUM EFFECTS AFTER DECOHERENCE 

We solve Eq.(2) numerically using a fourth-order spec- 
tral algorithm (numerical checks included carrying out 
simulations at different spatial and temporal resolutions). 
We chose A = O.f, D = 0.3 and set ^ = 1.0 initially. 
In order to understand the effects of the change in the 
mass term on the decoherence process we look first at the 
evolution of the quantum superposition of two Gaussian 
wave packets: 



where 

*i,2(z) = N(t)exp 



V{x,t = 0) = *i(ar) + V 2 (x) 
(x T L ) 2 



25 2 



(4) 

exp(±iP aO- (5) 



The initial W T consists of two Gaussian peaks W*' 2 sep- 
arated by a distance L (we chose L = 2.0 and Po = 0) 
and an interference term W,! nt . This quantum initial 
state has been widely used in the literature to illustrate 
decoherence phenomena (see [2] or [6] for example) and 
its evolution will make clear the physical nature of the 
effects we will observe. In the next Section we will chose 
a more realistic initial condition in terms of the dynamics 
of a phase transition. 

In order to visualise deviations from classicality effec- 
tively we define the auxiliary quantity [8] , 



T(t) 



dxdp [\Wr\-Wr]. 



(6) 



When the Wigner distribution is positive and possibly 
identifiable with a classical probability distribution, T is 
zero. However if T is positive W r must have negative 
values due to quantum interference terms. We can thus 
use positivity of the T function as sufficient condition for 
non-classical behaviour. 



A. Mass quench at constant temperature 

We start the simulation by evolving W r for some time 
with the positive mass squared potential, in the pres- 
ence of the bath. During this period, the initial quantum 
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interference terms are quickly damped by the environ- 
ment. Thus, for an early time tD 1 , the system decoheres 
and one is able to distinguish two classical probability 
distributions corresponding to the two initial Gaussian 
peaks evolving over phase space. Suddenly, at t — t c 
we change the frequency of the system from the initial 
positive value fig to a final — fig. The evolution picture 
changes dramatically when the frequency becomes neg- 
ative and instabilities arc introduced in the system. In 
Fig.l we can see the behaviour of T. Starting from a large 
initial value, T quickly tends to zero as quantum fluctu- 
ations vanish and the systems becomes classical. The 
potential is quenched at t c and shortly after the system 
displays once again quantum behaviour for a period of 
time. 
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FIG. 1. Evolution of T when the potential changes its fre- 
quency from Q,q — 1 ^ —1. D = 0.3 and A = 0.1. 

In order to understand this process we go back to early 
times, before the quench. From t — up to t = t c the 
diffusion coefficient D causes the system to decohere, de- 
stroying quantum interference terms in a time that can be 
estimated to be of the order of t-o 1 ~ 1/ (4LqD) , where Lo 
is the initial space separation between the peaks of the 
Gaussian wave packets (see [2]). The normal diffusion 
term is dominant with respect to the quantum correc- 
tions, and thereafter the evolution is given essentially by 
the classical Fokker-Plank flow. For our choice of ini- 
tial conditions we have tn 1 ~ 0.2. This is roughly the 
time quantum interference terms in the Wigncr function 
should fall to 1 je of their initial value (we have checked 
that this is compatible with the decay of T in the ini- 
tial period of evolution in our simulations). As soon as 
the frequency becomes negative, an unstable point forms 
in the centre of the phase space with associated stable 
and unstable directions. These are characterised by Lya- 
punov coefficients A with negative and positive real parts 
respectively [9]. 

The new type of dynamics gives rise to the possibility 
of squeezing along the stable direction. The exponential 
stretching of the Gaussian packets in one of the directions 



due to the hyperbolic point is compensated by an expo- 
nential squeezing. This will lead to a growth of gradients 
in the Wigner function that will make the quantum term 
in Eq.(2) comparable to the others. As a consequence 
the system will be forced to explore the quantum regime 
again. In a more quantitative fashion we have that the 
time dependence of the package width in the direction 
of the momenta after the quench is given by a p (t) = 
<7 p (t c ) cxp [— A(t — t c )], where a p (t c ) is the corresponding 
width at the time in which fio changes sign. From this 
we can estimate the p— derivatives of the Wigner function 
to grow as <9™W r oc a~ n (t c ) exp [nA(t — t c )]W r . Clearly 
higher order derivatives grow faster and at some point the 
quantum term with its third order derivative will be of 
comparable magnitude to the classical terms in the Pois- 
son brackets (which are first order). This will happen 
(see [9]) when the ratio dpW T /d p W T becomes of the or- 
der of x 2 = d x V(x)/9%V(x) ~ fig/A which characterises 
the scale of nonlinear terms. From this the time at which 
quantum effects become relevant is calculated to be 



t x <~ t c + A 1 In [x<J P {t c 



(7) 



In the simulation used in our example we chose t c — 2.5 
(later than the time when the Wigner function becomes 
definite positive). We evaluate x ~ 3.2 and numerically 
estimate a p (t c ) ~ 2.7. We are also assuming the Lya- 
punov coefficient is given by the value corresponding to 
a linear potential A = 2fig = 2.0. Therefore, the time 
in which quantum effects start being relevant is given by 
t x ~ 3.8. This is in good agreement with the time at 
which the Wigner function displays negative values once 
again, as can be seen in Fig.l. 

From this point onwards quantum contributions in- 
crease, their growth being limited by diffusion effects 
which limit the squeezing of the Wigner function. The 
bound on the width of the packs is given by cr c = ^/2D/A 
[2,9]. We use this to estimate the second decoherence 
time scale. We assume that quantum effects become 
maximal at a certain t max (when in the numerical simula- 
tion r reaches its maximum) with a corresponding pack 
width a p (t max ) and that decoherence is effective after the 
time when squeezing becomes of the order of the limiting 
value. This implies 



max 

+ A 1 ln[o- p (i max )/cr c ], 



(8) 



which defines the decoherence time after the critical time. 
Using CT p (i max ) ~ 4 and cr c = 0.5 we obtain £d 2 ~ 6.3, in 
reasonable agreement with the simulation time for which 
quantum effects arc exponentially suppressed (see Fig.l). 



B. Mass quench with changing temperature 

The pattern of classical-quantum-classical behaviour 
found in the above system with explicit time dependence 
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is observed in more generic situations. As a second ex- 
ample we have solved Eq.(2) allowing the bath tempera- 
ture to decrease simultaneously with the change in sign 
of the frequency term. These conditions take us some- 
how closer to what would happen in a true second-order 
phase transition caused by a temperature quench. As a 
consequence, the diffusion coefficient, proportional to T, 
goes at t c from an initial high temperature value Do up 
to a final lower value Df (still in the high temperature 
regime in order to ensure the validity of Eq.(2)). In Fig. 2 
we see the effect of changing the temperature with the 
classical potential (except for D all simulations param- 
eters are the same as in Fig.l). The analysis used in 
the previous example can be easily reproduced for this 
case. Both the initial decoherence time tu 1 and the time 
for the re-introduction of the quantum fluctuations t x 
remain unchanged as they do not depend on the temper- 
ature of the environment. The second decoherence time 
£d 2 is larger for a weaker diffusion term (we have used 
D f = 0.1 in Fig.2-a and D f = 0.003 in Fig.2-b). We 
have obtained respectively tu 2a ~ 6.5. and tu 2b ~ 7.4. 
In the lowest temperature case (Fig.2-b) the analytical 
prediction matches the numerical result poorly. 
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FIG. 2. Evolution of F when the potential changes its fre- 
quency from Qq = 1 — f — 1. A changing in the environment 
temperature is considered D — 0.3 — > 0.1 in curve (a) and 
D = 0.3 -> 0.003 in (b). 

This is due to the fact the estimation does not take into 
account the oscillations in the rate of decoherence com- 
ing from different orientations of the interference fringes 
when the Wigner function is moving around the unstable 
point. As the diffusion coefficient is smaller, the second 
decoherence time grows and the approximation of the 
upside-down potential in no longer valid. In any case, 
the analytic result can still be used as an estimated lower 
limit for the second decoherence time. We have included 
it in our analysis in order to emphasise how dramatic the 
quantum effects are during the quenched transition. 

It is helpful to look at the Wigner function directly 
in order to further clarify which regions of phase space 



are responsible for turning T positive. In Figs. 3 we show 
W(x,p, t) for the quench case corresponding to Fig.2-a. 
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FIG. 3. Stroboscopic phase space for the evolution in 
Fig.2-a. Horizontal axis corresponds to x, vertical axis to 
p. The medium grey shade on the background corresponds to 
zero values for the Wigner function, lighter and darker shades 
respectively to positive and negative values of W(x,p) 

The four plots in the left column correspond to the 
decoherence period before the quench. The two Gaus- 
sian peaks (light spots) rotate in phase space around the 
minimum of the potential while the negative components 
(dark patches) of the Wigner function are cleared away 
by the environment. When the potential changes (right 
column) the wave packets start spreading and exploring 
the new non-linear regions of phase space giving rise to 
the dark interference patches. For longer times decoher- 
ence takes over again and the Wigner function becomes 
once more positively defined. 
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IV. SINGLE INITIAL GAUSSIAN STATE 

As a further example we take a single Gaussian state 
centred at the global minimum of the quartic potential 
as initial condition. This is a more reasonable initial 
condition in terms of a realistic phase-transition, mim- 
icking a high-temperature thermal distribution. It will 
also allow us to see that the above results are not an ar- 
tifact of the initial state. This initial Wigner function 
is already classical and so we ignore the initial evolution 
period and take t c = 0. Fig. 4 and Fig. 5 show the T func- 
tion for the same quenches as before (without and with 
temperature change respectively). The initial classical 
configuration (r = 0. for the initial time) develops quan- 
tum effects as the classical potential and the temperature 
change. The relevant time scales are evaluated as before 
and once again, the estimates are in good agreement with 
the simulation results. In the constant temperature case 
(7 p (t c = 0) ~ 0.7 which gives t x ~ 0.6 (see Fig. 4). We also 
have <Jp(t max ) ~ 3.2. and a c ~ 0.5 leading to tu ~ 3.2, 
which agrees with the numerical result. 
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FIG. 4. Evolution of T when the potential changes its fre- 
quency from Qq — 1 — > — 1 for one Gaussian initially centred 
at x = 0. (A = 0.1). 

Fig. 5 shows the cases where the change in frequency is 
followed by a change in the environmental temperature 
(same coefficients as in the example of Fig. 2). For Fig.5-a 
Cp(imax) ~ 3.3 and a c = 0.3, and therefore the decohcr- 
ence time is tr> ~ 3.5. This scale is in good agreement 
with the numerical result. The estimation for Fig.5-b 
gives a decoherence time tu ~ 4.5 which again (as in the 
case of Fig.2-b) fails to fit the numerical result. We have 
included it in our analysis in order to emphasise how 
dramatic the quantum effects are during the quenched 
transition. 
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FIG. 5. Evolution of V when the potential changes its fre- 
quency from Qq = 1 — * — 1 for one Gaussian initially centred 
at x — 0. (A = 0.1). A changing in the environment tempera- 
ture is considered D — 0.3 — > 0.1 in (a) and D — 0.3 — ► 0.003 
in (b). 



V. LINEAR ENTROPY 

One of the most salient features of the quantum to 
classical transition concerns the production of entropy 
as a consequence of the entangling interactions between 
the system and the environment. In order to clarify the 
nature of the post-decoherence quantum effects in the 
systems simulated above we have looked at the corre- 
sponding time evolution of the linear entropy. This is 
given in terms of the density matrix by [4] : 



S ; (t) = -ln[Tr[p 2 (t)]]. 



(9) 



This quantity can be easily obtained from the Wigner 
function giving a good measurement of the 'loss of purity' 
of the system as it interacts with the bath (see [2,4]). We 
found that as expected the entropy increases through- 
out the whole evolution. The system starts as a pure 
state and while interacting with the heat bath it looses 
coherence and simultaneously starts behaving as a clas- 
sical ensemble. When the potential changes it evolves for 
some time as a quantum mixed system but the original 
'purity' is never recovered. In this sense the decoherence 
process is irreversible. In terms of the Wigner function 
the linear entropy is related to the area of its non-zero 
component in phase space. Due to the coupling to the 
environment the total area is not conserved, the Wigner 
function keeps spreading at all times leading to perma- 
nent growth of the entropy. 

Fig. 6 shows the time dependent linear entropy (top 
plot) and its production rate (bottom plot) for the two 
different initial conditions considered before, in a quench 
with fixed environment temperature. 
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rate. 

We should stress that a growing linear entropy func- 
tion does not imply classicality (positivity of the Wigner 
function is an extra necessary condition in order to have 
a classical probability distribution). Increase of Si tells 
us that the pure initial quantum state is evolving into a 
mixed state. It does not of course, tell us whether this 
mixed state is a classical or quantum one. In particular 
this is the case between the quench time and the second 
decoherence time. During this period quantum effects 
are re-introduced while the linear entropy is still growing 
(faster even since its production rate increases). 

VI. FINAL REMARKS 
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FIG. 6. In the top figure, we show the linear entropy for 
the de-localised (solid line) and single Gaussian (dashed line) 
initial states. In the bottom figure, we have the entropy pro- 
duction rate vs time for the same initial conditions (solid and 
dashed lines as above). 

In the case of the double-gaussian initial state (solid 
line) there is an initial period of evolution up to the first 
decoherence time tD x = 0.2, where the linear entropy 
grows as a consequence of diffusion effects (as during 
the whole evolution) and also due to the disapearence 
of initial interference terms which are washed away by 
the environment. As these vanish the entropy produc- 
tion rate decreases as can be seen in the bottom plot. 
After an oscillation in the rate caused by the rotation of 
the Wigner function in the phase space which generates 
some low amplitude interference terms, the rate reaches 
a minimal value near the quench time at t = t c . After t c , 
the entropy rate starts growing again as the system gets 
rid of the newly induced interference terms. Finally at 
i_D 2 , the entropy rate decreases to a low, slow decaying 
value driven by diffusion only. 

The single gaussian evolution (dashed line) confirms 
this picture. From a low initial value (the initial state 
is free from negative terms) the entropy production rate 
grows as the quench generates interferences. Later, the 
environment cleans them out leading to the final decaying 



We have shown, using an exact numerical evaluation 
of the Wigner function that quantum effects can be re- 
introduced after decoherence in several systems with ex- 
plicit time dependence. These quantum effects are origi- 
nated when the changing dynamics introduce instabilities 
in previously stable regions of the phase space. When this 
happens the dynamics of the Wigner function becomes 
more relevant than the decoherence effects due to the 
environment (and the lowest the final bath temperature 
the more dominant these are). The system then displays 
quantum behaviour for a length of time until the envi- 
ronment manages to catch up and force classicalization 
once again. 

Since all examples so far were based on systems de- 
scribed by a double well potential one could wonder 
whether our results could be a consequence of possible 
tunnelling phenomena between the two minima. Tun- 
nelling is possible between symmetry related eigenstates 
with energy below the barrier. The tunnelling time-scale 
for each pair is well known to be inversely proportional 
to the energy splitting of the symmetry related pair of 
eigenstates. For the parameters of our system (ft 2 = 1, 
A = 0.1, h = 1) only seven pairs of states are found 
below the barrier. Their energy splittings range from 
AE ~ 1CT 12 to AE 7 ~ 10~ 2 , and thus the tunnelling 
would firstly be expected after t ~ 100. Therefore and 
considering the time scales in which our simulations take 
place, tunnelling should play no role. The stretching and 
folding of the Wigner function responsible for the ob- 
served effects happens on both 'sides' of the potential 
well independently. This is in agreement with the con- 
clusion invariably found in the literature (see for example 
Ref. [10]) that tunnelling takes place rather slowly when 
compared with all natural time-scales in the system. 

In order to confirm directly that tunnelling phenom- 
ena are not responsible for the effects observed, we solved 
numerically the problem of a single Gaussian packet can- 
tered at x — 2 evolving in the usual quartic potential 
but with its motion restricted to x > 0. The resulting 
T is shown in Fig. 7. As before quantum behaviour is 
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swiftly recovered, the corresponding time scales being in 
good agreement with the analytical estimates. 
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FIG. 7. Evolution for a single Gaussian packet cantered at 
xo = 2 (t — 0) in a x > single-welled potential (and an 
infinite wall at x = 0). We found t x = 0.8 and to = 7.4 
(fig = -1 and D = 0.01). 

Our results open up several interesting possibilities. 
The most obvious one would be to try to 'maximise' the 
recovering of quantum effects to the extent of making 
them effectively permanent. An oscillatory frequency [11] 
that would continuously force instabilities into the sys- 
tem could prevent classicalization or at least postpone it 
for a great length of time. 

In terms of the specific case of the dynamics of a 
second-order phase transition one could expect quantum 
effects to be present. Though the model used is a crude 
simplification of what happens in a realistic phase tran- 
sition, the same features of time-dependent introduction 
of non-linearities would be present in that case, leading 
to similar, probably stronger quantum effects. Critical 
properties of infinite dimensional systems such as crit- 
ical slowing down could play an interesting role in the 
process. 
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